Complete Set of Homogeneous Isotropic Analytic Solutions 
in Scalar-Tensor Cosmology with Radiation and Curvature 

Itzhak Bars*, Shih-Hung Chen^, Paul J. Steinhardt* and Neil Turok^ 

* Department of Physics, University of Southern Cahfornia, Los Angeles, CA 90089, USA 
^Perimeter Institute for Theoretical Physics, Waterloo, ON N2L 2Y5, Canada 
"^Department of Physics and Princeton Center for Theoretical Physics 
^ ! Princeton University, Princeton, NJ 08544, USA 

2 ■ Abstract 



^3 

a^: 



X: 



We study a model of a scalar field minimally coupled to gravity, with a specific potential 
energy for the scalar field, and include curvature and radiation as two additional parameters. 
Our goal is to obtain analytically the complete set of configurations of a homogeneous and 



p [ isotropic universe as a function of time. This leads to a geodesically complete description of 

^ ■ the universe, including the passage through the cosmological singularities, at the classical 

level. We give all the solutions analytically without any restrictions on the parameter space 
^ ■ of the model or initial values of the fields. We find that for generic solutions the universe 

^ . goes through a singular (zero-size) bounce by entering a period of antigravity at each big 

^ _ crunch and exiting from it at the following big bang. This happens cyclically again and 

^ ■ again without violating the null energy condition. There is a special subset of geodesically 

complete non-generic solutions which perform zero-size bounces without ever entering the 
antigravity regime in all cycles. For these, initial values of the fields are synchronized and 
quantized but the parameters of the model are not restricted. There is also a subset of 
spatial curvature- induced solutions that have finite-size bounces in the gravity regime and 
never enter the antigravity phase. These exist only within a small continuous domain of 
parameter space without fine tuning initial conditions. To obtain these results, we identified 
25 regions of a 6-parameter space in which the complete set of analytic solutions are explicitly 
obtained. 

PACS: 98.80.-k, 98.80.Cq, 04.50.-h. 
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I. INTRODUCTION 



Scalar-tensor theory has been one of the most popular tools for building models in cosmol- 
ogy. It is sufficiently simple while havinga variety of applications. One common application 
is to describe early universe inflation [l-3| where the scalar plays a central role in driving a 
period of accelerated expansion that solves the homogeneity (horizon) and the flatness prob- 
lems, and generates the primordial density perturbation that seeds the subsequent growth 
of the large scale structure ^-l3|. Alternatively, some ekpyrotic models for the early universe 
models 8|, |9|] also utilize scalar-tensor theories to produce a period of a slow contracting 
phase before a big crunch that eliminates the horizon problem, and solves the flatness prob- 
lem as well. With some speciflc matching rules inspired by a colliding two branes picture, 
ekpyrotic models can also generate scale invariant density perturbations as observed today. 
Another application of scalar-tensor theories is to produce the late time acceleration of the 
universe that is inferred from type lA supernova observations [lo|. Such models are called 
"quintessence" models where a small nonzero dynamical vacuum value of a scalar potential 
replaces the cosmological constant In addition, by using a conformal transformation, 
it has been shown that modifled gravity theories, such as / (R) gravity, are equivalent to 
scalar-tensor theory with a speciflc scalar potential 12 1. 

Despite such broad applications of scalar-tensor theories, only a few isolated examples 
of analytic solutions have been found so far. This is because the coupled second order 
nonlinear differential equations are hard to solve analytically. Our goal in this paper is to 
provide a full set of analytic solutions that give all possible conflgurations of a homogeneous 
and isotropic universe as a function of time. This expands on our previous work in 13l-ll5| 



by including additional degrees of freedom, in particular radiation. The effects of anisotropy 



are discussed elsewhere 
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17|. 



Our overall approach in this paper is in contrast to speciflc analytic, approximate or 
numerical solutions that are usually flne tuned from the point of view of initial conditions 
and/or the potential energy function V (a) for the scalar field, to force a solution in which 
the universe has a particular desired behavior as motivated by prejudices and observations. 
Instead, we would like to understand the global structure of solution space that can emerge 
from a class of theories, so that we can gain a better understanding of how the features 
of our own universe could emerge. Obtaining the full set of classical solutions can provide 
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some such insights. Indeed through our solutions we gain new understanding about general 
generic behavior as we will see below. 

We can obtain the full set of analytic solutions of the scalar-tensor theory for several forms 
of the potential energy function V (a) for the scalar field. In this paper we concentrate on 
the case: 



6^^ 



csinh^(^^-aj+6cosh^(^^-ajJ, (1) 

where the parameters b and c are dimensionless free parameters, and is the reduced 
Planck mass = sj^^ = 2.43 x 10^^^^. We note that this potential has familiar 
features. For example, if {h + c) > 0, depending on the various values and signs of b, c, 
V{a) has a single well or double well with stable minima, similar to other potentials used in 
cosmological applications. Because our analysis has a broad range of applications beyond 
cosmology, we will classify all the solutions regardless of the physical application. In other 
papers, including [l6, the role of these solutions is discussed in a cosmological setting. 
Some of these cases will be pointed out briefly later in the paper. In a separate paper we 
will present the analytic solutions for the potentials V (a) = ^he^v'^'^/^ for arbitrary p and 
V (cr) = ^ ^^g"2ft(T/v^ _j_ gg-4Kcr/v^j ^ where 6, c, p are dimensionless real parameters. 

For the potential [H we can solve the Friedmann equations exactly for all time intervals 
before or after the big bang. The method, which is based on conformal symmetry, was 

n n 

introduced in previous papers [13[-[17[. It was applied to the cases of the flat and curved 
isotropic Friedmann universes without radiation or matter 1^. It was also applied to the 
case of an anisotropic universe in the absence of the potential energy, but with the inclusion 
of radiation [17]. In this paper, we further generalize our method for the isotropic universe to 
include both curvature and radiation with the potential, where radiation is taken in the form 
of a perfect fluid. The inclusion of radiation is a simple mathematical exercise beyond our 
previous paper 1^, but it describes richer physics, and leads to more complicated analytic 
expressions for the solutions. So, the reader may wish to flrst understand the previous work 
in a simpler setting 

Including radiation, the model is deflned by 4 parameters, namely (6, c) in the potential, 
the curvature parameter K in the metric, and pr-, the energy density of radiation when the 
scale factor is a = 1. The two flelds of interest are the time dependent scale factor a (r) 
and the scalar fleld cr (r) . Their initial conditions (a (tq) , d (tq) , a (tq) , a (tq)) introduce 4 
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more parameters which enter the general solution to second order. However, the zero energy 
condition including gravity (or Friedmann equation) eliminates one of the initial values, and 
one other initial value can be absorbed into a redefinition of the initial time tq by using the 
time translation symmetry of the differential equations. Hence the complete set of solutions 
are described by 4+2=6 parameters given by {b,c, K, pr, E, a (tq)) , where a single energy 
parameter E is used conveniently instead of the two related initial velocities a (tq) , a (tq). 
We will give all the analytic solutions without putting any conditions on these 6 parameters 
at any tq. It turns out that there are 25 distinct regions of this parameter space in which 
the solutions take different analytic forms in terms of Jacobi functions. The 25 regions and 
the corresponding solutions are given explicitly in the Appendix. If an unstable potential 
with {b + c) < is of interest there would be more regions and corresponding solutions; with 
appropriate modifications these can be obtained from those available in the Appendix. 



As in our previous work Il5| . we find that in the Einstein frame the generic solutions 
are geodesically incomplete at the cosmological singularity (see Eqs. (3) and (4) in jl^). 
We will then find that the knowledge of the full set of solutions suggests how to complete 
the space and make it geodesically complete for the generic solution. This completion 
includes time intervals during which the gravitational field effectively acts like a repulsive 
force (antigravity) rather than an attractive force (gravity). So, the generic solution has 
alternating time intervals of gravity and antigravity as graphically illustrated in Fig. [2l In 
this figure, each time the trajectory crosses the 45° lines, the universe transits from gravity 
to antigravity or vice-versa. There is however a subset of non-generic solutions that are 
geodesically complete in the Einstein frame without ever entering the antigravity region. 
These are of two types: (i) singular zero-size bounces at the cosmic singularity without 
violating the null energy condition as shown in Fig.|3]and (ii) non-singular finite-size bounces 
as shown in Figs. (23,24) in jl^. 




The zero-size bounces (i), which are classified in tables I, lib and HI, with the related an- 
aljd;ic expressions in the Appendix, are obtained by synchronizing and/or quantizing some 
initial values. These tables provide the most general parameter subspace (within the 6- 
parameter set) for which the geodesically complete singular bounce occurs without anti- 
gravity. The parameter subset consists of 4 continuous and one quantized parameter, which 
is obviously smaller than the available continuous 6-parameter set. Despite the fact that 
this subset of solutions may be considered a set of measure zero as compared to the full set. 
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it is distinguished as the only zero-size bounce set of solutions that are geodesically complete 
in the Einstein frame, and do not enter the antigravity region at any time in any cycle. 

The finite-size-bounces (ii), which are classified in table Ila, with the related analytic 
expressions in the Appendix, describe a universe that contracts up to a minimum non- 
zero size, at which point the spatial curvature causes the universe to bounce back into 
an expanding phase. This kind of spatial curvature-induced bounce is already familiar in 
cosmology. Here we provide analytic solutions for the finite-size bounces. As the universe 
turns around to repeat such cycles, the minimum size is not necessarily the same in each 
cycle, as this depends on the parameters. Such solutions occur when the parameters satisfy, 
Pr < K^/lQb and 0Ln (^o) > K/Ab > s^^^ (tq) (see Table Ila). Note that there are still 6 
parameters, so this is a continuous set, but it is a restricted region of parameter space or 
initial values. 

The solutions above - the generic case, type (i) or type (ii) bounces - are the exact and 
complete set of solutions in the absence of anisotropy. Although anisotropy can be neglected 
as the universe expands, it can grow to be a dominant effect near the singularity. We do 
not consider those cases here; they are described in la, ll7|] where it is proven that there 
is an attractor mechanism that is independent of initial conditions. The attractor distorts 
the zero-bounce solutions, both the generic or non-generic type (i), to behave in a unique 
way such that these solutions must undergo a big crunch/big bang transition by contracting 
to zero size, passing through a brief antigravity phase, shrinking to zero size again, and 
re-emerging into an expanding normal gravity phase. 

This paper is organized as follows. In Sec. [Tll we introduce the standard scalar-tensor 
theory, with a single minimally coupled scalar field a (x) , as the gauge fixed version of a 
locally scale (Weyl) invariant reformulation of Einstein's theory of gravity that contains two 
conformally coupled scalar fields (x) , s (x) . This "Weyl lifted" version has an extra scalar 
field as well as a local Weyl symmetry that compensate each other, so that the physical 
degrees of freedom are the same number as in the standard formulation of the theory. This 
model and method of solution emerged directly from the 2T-physics formulation of gravity 
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20|. The model was also inspired in the context of braneworld notions 2l|, |22|] that led 



to the colliding brane scenario for cosmology [8|, |9|, |23|, |2J|. Recently 'tHooft also motivated 
the same Weyl invariant theory because of its ability to give a better description of black 
and white holes in a convenient gauge 3]. Such gauge choices, including the E,'y,c and 
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s gauges discussed in this paper, are just a small subset of examples of 3+1 dimensional 
"shadows" that 2T-physics yields as dual forms of the same parent theory that unifies them 
in 4+2 dimensions. 

The Weyl lifted version has no dimensionful constants, not even the gravitational con- 
stant. The extra scalar field can be eliminated by gauge fixing it to a dimensionful constant, 
thus reaching what we call the "Einstein gauge" (ii^-gauge), the standard formulation of the 
theory in the Einstein frame with the usual gravitational constant and the scalar field a. 
The advantage of the Weyl lifted version is that it allows us to choose another more conve- 
nient gauge that we call the 7-gauge, in which the cosmological equations greatly simplify 
and can be solved analytically. The full set of solutions is then mapped to the Einstein 
frame by a gauge transformation from the 7-gauge to the i?-gauge, and verified that they 
are the solutions of the Friedmann equations. In this process we only add gauge degrees of 
freedom to Einstein's theory. But in the presence of the gauge degrees of freedom we find 
naturally the geodesically complete space and understand much more clearly the nature of 
the complete space of solutions. In particular we learn that geodesic completeness requires 
an extension of the domain of the original fields in the scalar-tensor theory in the Einstein 
frame, such that with this extension, the same fields can describe also an antigravity region 
not captured at first sight. 

In Sec. Illll we show how the complete set of solutions of the Friedmann equations, includ- 
ing curvature and radiation, are obtained analytically without constraining the 6-parameter 
space. The complete set of solutions is explicitly given in the Appendix, where in 25 different 
regions of the parameter space the analytic expressions take different forms. These solutions 
are all cyclic and geodesically complete in an enlarged domain as described above. 

The non-generic solutions of type (i) with zero-size bounces, and type (ii) with finite-size 
bounces, which never enter the antigravity regime, are still geodesically complete in the 
gravity domain. In Sec. IIVI we determine the constraints on parameter space and initial 
conditions that distinguish the geodesically complete solutions in the restricted Einstein 
frame. The corresponding parameter spaces are classified in Tables I, IIa,b and III. 

In Sec. |V] we comment on the generic solutions that are geodesically complete provided 
an anti-gravity regime is included. 

In Sec. I VI I we summarize our results. In an Appendix we list all the analytic cosmo- 
logical solutions predicted by our model introduced in Sec. HIl 
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II. THE WEYL LIFTED MODEL 



Our approach begins with the standard action typically used in cosmological models that 
describe a scalar field a {x^) minimally coupled to gravity 



i^R {qe) - \gE9^i(ydua - V {a) 



^radiation + matter 

We are able to solve for the complete set of homogeneous, isotropic, cosmological solutions 
of this model when the potential V {a) is given as in Eq. ([1]). Through these solutions we 
discover the geodesic completion of the space through the cosmological singularity. 

To solve the equations and to understand the geodesic completion we will use a device that 
we call "Weyl extension" in which the model is enlarged by adding gauge degrees of freedom 
that are compensated with a local scaling (Weyl) symmetry. The local scaling symmetry 
does not allow the usual Einstein- Hilbert term R^Qe) /2k^, but allows conformally coupled 
scalars. The following action, which will be shown to be related to the one above by a gauge 
choice, contains two conformally coupled scalars, 0, s, interacting with the curvature term 
with the coefficient dictated by the gauge symmetry 

S = j d'x^ [lg>^'^d,<pd.<p - ^g^^d.sd^s + ^ (0' - s') R {g) - <p'f ) . (3) 

The function f{z) is determined by the scalar field potential (see below) and, in general, 
can be an arbitrary function of the gauge invariant ratio ^ = | and still maintain the Weyl 
symmetry. The metric g^y can differ from the metric g^^ in Eq. (|2]) by an arbitrary conformal 
factor because of the gauge symmetry under the following local transformation 

(j) -)■ VLcj), s -)■ VLs, g^u Vr'^g^y, 

where VL [x) is an arbitrary function of spacetime. 



This action was described as the conformal shadow of the 2T-Gravity action [19|, l20| in 
4+2 dimensions. In this setting, it was shown that the local scaling gauge symmetry is a 
remnant of general coordinate transformations in the extra 1+1 dimensions. 

We draw attention to the fact that has the wrong sign kinetic term while s has the 
correct sign. It appears as if is a ghost; however, since the ghost disappears for some 
choices of gauge (e.g. the gauge that restores the Einstein gravity form of the theory), there 
is no real problem with ghosts or unitarity. In this connection, note also that there is no 



gravitational constant; rather, the factor ^ (0^ — s^) effectively behaves like a gravitational 
parameter which replaces the usual expression (leTrG)"^ = (2k2) where G is the Newton 
constant. If (j) had the opposite sign kinetic term, then this factor would become purely 
negative (—0^ — s^) in order to maintain the Weyl symmetry, but then the gravitational 
parameter would have the wrong sign. This is the reason why cf) must be introduced initially 
with the "wrong" sign, so that the gravitational parameter ^ (0^ ~ ■5^) is positive at least 
in some regions of field space. 

Fermionic and gauge fields, as well as more conformally coupled scalars, can be added, to 
construct a completely Weyl invariant model, such as the Standard Model of particle physics. 
The Higgs field must also couple as a conformal scalar to preserve the Weyl symmetry. 
The Higgs mass term is not allowed by the gauge symmetry, but it can be generated by 
coupling the Higgs doublet H to the singlets 0, s in Weyl invariant quartic terms of the form 
{acl? + There are various possible models for the effective Higgs mass by choosing 
the parameters a, /3. One possibility is related to the fact that (0|; — s\) gets fixed to a 
constant in the Einstein gauge of Eq. (I6l). Another possibility emerges from our solutions, 
by noting that cosmologically Se = ^ sinh {^^^^ evolves to a field Se ~ cr (r) which is 
much smaller than the Planck scale, and at late times behaves almost like a constant that 
ca. the H,ggs ,„ass te™. This „.gM a.so explain the „,ass h,e..a..chy HQ- 

We discuss here four gauge choices that are useful for finding solutions and interpreting 
them: Einstein gauge (ii^-gauge), the 7-gauge, the supergravity gauge (c-gauge) and the 
string gauge (s-gauge). To distinguish the fields in each gauge we denote them by the 
subscripts E, 7, c, s respectively. We also define the following Weyl gauge invariant quantity 
X which plays an important role in our discussion 



X^^{-9)k<P'-s'). (4) 
Another gauge invariant is the ratio 

E-gauge: The usual Einstein gravity in Eq. ([2]) is obtained in the i?-gauge, in which we 
denote the fields by 0^; (x) , s^; (x) , (x) with a subscript E. In this gauge the gravitational 
parameter is gauge fixed to the usual Newton constant for all spacetime 

^^{<^l{x)-sl{x)) = {2^y\ (5) 
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Then parameterizing this gauge with a single scalar field a, 

(j)E {x) = ±^ cosh y-^^j ' (x) = ±^ sinh [-^J^j ' (6) 

we find that the conformally gauge invariant action yields the familiar action in Eq. ([2]) 
with an Einstein-Hilbert term {qe) and a minimally coupled scalar field cr, just as in 
Eq. ([2]). The gauge invariant x iii the i?-gauge becomes 
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X = {-9e)-^- (7) 

In a cosmological solution with the metric 

dsl = a| (r) {-dr^ + dsf) (8) 

where the 3-dimensional metric ds^ will be discussed later, {—gE)^ is just the scale factor 
of the universe, x = (^"e ('^) 5 so x niust vanish at the cosmological singularity at the time 
of the big bang ge {jb) = 0. This shows that the big bang corresponds to x (Jb) = 0, and 
since x is gauge invariant, the cosmological singularity is at x {'^b) = in all gauges. 

Note that the geometry completely collapses at the singularity at time tb in the £'-gauge 
since {—Qe (tb)) = 0, and this is why geodesies are incomplete. However, we will see that 
the geometry does not collapse at all at x (Jb) = in the other gauges, and this is how we 
are able to complete the geodesies and the geometry. 

'j-gauge: In the 7-gauge we choose (— 5'-y)3=constant for all r. Since the cosmological 
FRW metric in the E'-gauge is conformally fiat (even when the curvature is non-zero), 
ds"^ = a\ (r) (— e^dr^ — ds\p^yY)i we can at first discuss the general conformally flat metric 
which can always be put into the form ds^ = {x^) rj^y, where rj^^i, is the Minkowski metric. 
A conformally fiat metric becomes fully fiat in the 7-gauge ds"^ = rj^,^, namely the gauge 
choice (x) = 1, leading to R^g-y) = 0. Hence the degrees of freedom in the 7-gauge are 
0^, in a fiat Minkowski geometry. The gauge invariant x ^low takes the form 

xixn = ^{^',-s'^){xn, (9) 

and the gauge invariant action in ([3]), taken with any conformally fiat metric, reduces to a 
gauge fixed action in fiat Minkowski space 

'1 



(10) 



where (x) is a ghost since it has the wrong-sign kinetic energy. However, this ghost is 
removed as follows: before choosing the 7-gauge we recall that there was an Einstein equation 
for the metric, G^i, (g) = T^y, which must be imposed for any metric in the 7-gauge as well. 
For the conformally flat metric in the present case, the curvature vanishes in the 7-gauge, 
R^ivXa iv) = 0, leading to G^,^ {rj) = 0, and therefore T^^, (0^, s^) = 0. The vanishing of the 
stress tensor for the action is a constraint that removes the ghost 4>j. 

To include the effect of both anisotropy and curvature, the 3-dimensional part of the FRW 
metric, {dsDpj^^^, must be replaced by the corresponding 3-dimensional parts of the Kasner 
metric, Bianchi IX metric, Bianchi VIII metric, when K — 0, K > 0, K < respectively. 
We parameterize the Kasner metric {K = 0) as follows 

where ai {x'^) and a2 (x^) are the anisotropy fields that are taken to depend only on time in 
the homogeneous case considered in this paper. The Bianchi IX and VIII metrics contain the 
same ai 2 {x'^) , as well as the non-zero curvature parameter K that is included by generalizing 
{dx,dy,dz) to {dax,day,daz)- Both of these Bianchi metrics reduce to the Kasner metric 
when the curvature K vanishes. 

In the homogeneous limit, with an anisotropic cosmological metric of the form given 
above, ds"^ — €? (r) {—e^dr^ — dsg) , the Weyl invariant action reduces to the following effec- 
tive action for the homogeneous cosmological degrees of freedom 



5'efr= / dT 



(11) 



__L (a, (a</,))2 + X (a, ^as)f + i ((/'^ - .2) {a\ + di) 
-e aVV (f) +Pr - I (0^ - s^)a^^^(Q;i,Q;2) , 

where is the radiation density when the scale factor a = 1. The effective action is invariant 
under time dependent Weyl transformations 

a (r) ^ (r) a (r) , (0 (r) , s (r)) ^ ^ (r) (0 (r) , s (r)) (12) 

while ai 2 and e are Weyl invariant. Here v (cti, a-}) is the anisotropy potential which emerges 
from the curvature term (0^ — s^) R {g) 



V [ai, a2) 



1-4 sign {K) 



-4 sign {K) e-\/273««i cosh (\/2fi:a2) 



(13) 
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In the isotropic limit the anisotropy potential reduces to a constant linioi 2-s>o (^i, ^2) = K. 
For the Kasner metric the potential energy term is absent even if anisotropy is present since 
K = {}. 

For the cosmological solutions discussed in this paper we will concentrate on the isotropic 



case ai_2 — 0, and therefore v (ai, 0:2) — )■ K. 
geneous and anisotropic universes is given in 



he generalization of our discussion to homo- 



16| and [17|. 



In this action e (r) is related to the lapse function which is a part of the metric g^y (x) . 
Its presence insures r-reparameterization symmetry. The equation of motion with respect 
to e imposes a constraint on the degrees of freedom. This constraint is equivalent to the 
Goo = ^00 Einstein equation. One may choose a gauge for r, such that e (r) = 1, in which 
case the action simplifies. 

In the Einstein gauge of Eq. ([6]) the cosmological action is 

^-^ = I r 4 T. / A 6 2/ V ( ^ > 

-e [a%V [a) + p,, - ^^aj^v (ai, 02) 



In this action note that the qe (t) degree of freedom has the wrong sign kinetic energy term, 
and therefore it is potentially a ghost. However, the Hamiltonian derived from this action 
in the e (r) = 1 gauge is required to vanish as a result of the constraint, and this is just 
sufficient to compensate for the ghost. This constraint equation, which is called the zero 
energy condition, is the same as the first Friedmann equation. 
The 7-gauge is defined by fixing 

(r) = 1 (15) 
for all r. In the 7-gauge, the cosmological action is 



S2^ = I dr 



1 



(16) 



+ + S (0^ - 4) (al + al) 
-e [0^/ (s^/ ^7) + Pr - 1(07 - 4) ^ ("1' "2)] 
When V (cr) is given as in Eq. ([1]), then 0^/ (s^/0^) = 60^ + cs^ is purely quartic. In this 



action the 0^ degree of freedom has the wrong sign kinetic energy term. However, as in 
the case of in the previous paragraph, due to the r-reparameterization symmetry, and 
the corresponding zero energy constraint that follows from the equation for e, this potential 
ghost is ehminated. 

We can relate the -gauge and 7-gauge degrees of freedom to each other by a gauge 
transformation. More easily, we identify the gauge invariants x and s/0 in both gauges, and 
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from them we extract 

.2 



2 I 1 2 2 1 , 

^E = ^\<P,-s,\.a = —\n 



^^^^ 

From this, we see that, what appeared as a cosmo logical singularity in the i?-gauge at 
a\ (tb) = does not show up at all as a geometrical singularity in the 7-gauge since 
R^ivXu iv) = 0- Of course, the big bang of the i?-gauge must reflect itself again as the gauge 
invariant x (^b) = 0, which becomes (0^ — s^) (r^) — i- 0, however this is not a singularity of 
the equations within the 7-gauge. 

The absence of a geometrical singularity in the 7-gauge, and the simplicity of the equa- 
tions of motions for 0^, is the key for being able to solve all the cosmo logical equations and 
finding the complete set of solutions analytically. This is also what permits us to understand 
geodesic completeness. The geodesically complete geometry includes space-time regions of 
antigravity. In the antigravity regions the gravitational parameter ^ (0^ — s^) that appears 
in the general action becomes negative. The switching of the sign in generic solutions of the 
theory occurs precisely at x {'^b) = 0, which appears as the cosmological singularity in the 
i?-gauge, but this is a completely smooth point for the geometry in the 7-gauge, as well as 
all other gauges, except the £'-gauge. 

Supergravity gauge (c-gauge): In the c-gauge 13], we fix 0c = 0o where 0o is a constant for 
all x^. Then in ([3]) or (ITT]) there is only one scalar Sc (x) while the curvature term takes the 
form ^ (0Q — si (x)) R {qc (x)) . We see that the 0q term plays the role of the Hilbert-Einstein 
term while the overall structure ^ (0q — si) R (qc) is similar to that found in supergravity, 
including the Kahler potential. In fact, this model can be regarded as a toy model for the 
scalar sector of a full supergravity model. Then the term {—slR{gc)) allows us to identify 
si {x) as the analog of the Kahler potential in supergravity. The gauge invariant becomes 
X = ■^(— 5'c)^(0o ~ ^1 i^))^ ^-iid for a cosmological solution it takes the form 

X = y(0o-ScW)«cM- (18) 

Since we have all the solutions, we can verify that the cosmological singularity x {'^b) = 0, 
which in the Einstein frame is at as {tb) = 0, occurs when si {tb) = 0o, rather than a1 (r^) = 
0. Hence the metric gft" in this gauge is not singular at the Big Bang. This is because the 
quantity Sc/0o = 5^/0^ = S£;/0£; is gauge invariant and takes the value (s/0) (r^) — )■ 1 at the 
singularity in any gauge. (The reason that the i?-gauge ^ {(f)% (x) — s% (x)) = {2k^) ^ can 
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remain finite even at the singularity x (Jb) = 0, is because in the expression — s^) (r^) = 
01; (tb) (1 — s\/(t)\) {tb) , when the second factor vanishes, the factor (p^ {tb) blows up so 
that — s\) remains unchanged at all x^, including at the singularity at time tb-) This 
shows that, similar to ^ (0q — s^) R [Qc) , the effective gravitational constant in the curvature 
term in the action of a supergravity theory is typically expected to switch sign at the big 
crunch or big bang. So the phenomena we discuss here, including antigravity regimes related 
to geodesic completeness, are expected generically in typical supergravity theories. 
String gauge ( s- gauge): We also discuss the string gauge (s-gauge) to connect to the typical 
structures in string theory. The string frame in d dimensions is defined by the following 
form of Lagrangian 

^string gauge = ^2^'^" ids) + Q^," d ,^d,^ - ^ ($))^ , (19) 

Note the wrong sign and unusual normalization of the kinetic term for the "dilaton" $ . 
When the transformation from the string-frame to the E-frame is performed by the substi- 
tution 

{9s),. = ev^* i9E),, , $ = v^cr, (20) 
then the right sign and normalization of the Einstein-frame a field emerges 



-^Einstein gauge 



[^R {9e) - \9'Ed,od.<y - V (^)) . (21) 



The string gauge in Eq. (fT9|) . for (i = 4, is obtained from our Weyl invariant action (j3]) by 
choosing the following gauge in which 0s, Ss are expressed in terms of a single scalar $ (this 
is analogous to 4>e,se given in Eq.([n])) 

(ps (x) = — e"^*^"^^ cosh 

Ss (x) = ^^e~^^^^^ sinh 

while the metric g'^'^ is labelled with the letter s to distinguish it from the metric in another 
Weyl gauge. This shows how the dilaton field $ of the non-Weyl invariant string theory can 
be related to the fields 0, s in the Weyl invariant theory. 

Then the effective cosmological action in the string gauge takes the form of Eq.f llip where 
the gauge fixed form for (0^, Sg) is inserted. The remaining degrees of freedom in the string 
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gauge are then (os, a) ,where a is the same field as the one in the Einstein gauge except for 
the overall normalization in Eq.f l20|) . while is related to a^; simply by the transformation 
in Eq.(I2U]) for d = 4 

V2^a. (22) 



Therefore, by using the relation between the -gauge and the 7- gauge in ( ITTI) . we can relate 
the s-gauge degrees of freedom, a^, $, to the 7-gauge degrees of freedom as follows 

- 4) 



K 

~6 



(23) 



$ = 73 In 



The expressions given above are consistent with the gauge invariants X; (s/^) , a0, as, as 
expressed in the s, 7, E, c-gauges, such as a^^s = aE4>E = «c0c = ^707, etc. 

In solving the cosmological equations various gauge choices turn out to be useful. In 
particular, for the cases we have solved, the 7-gauge turned out to be the most useful. By 
using the relations between gauges displayed above, the solutions for (r) , (r) imply the 
solutions of the degrees of freedom in the other gauge choices. 



III. SOLVING THE EQUATIONS WITH RADIATION AND CURVATURE 

The crucial step introduced in [l3^, as suggested by 2T-physics, is to take advantage of 
the gauge symmetry of the conformally invariant action ([3]). The Einstein equations derived 
from the action ([3]), taken with only time dependent isotropic fields, yields gauge covariant 
cosmological equations for three fields a (r) , (r) , s (r) in any gauge. If one chooses the 
Einstein gauge given in Eq. ([6]) then the corresponding gauge fixed fields a^; (r) , cr (r) satisfy 
the Friedmann equations including radiation pr, as given below in Eqs. (12511271) . Instead, if 
one chooses the 7-gauge defined by (x) = 1, the remaining fields 0^ (r) , (r) turn out to 
satisfy the decoupled equations (I28ti30p that can be solved exactly. Then, by performing a 
gauge transformation that relates the two gauge fixed configurations we obtain the full set 
of solutions of the Friedmann equations. 

The field re-definition from a, to 07, S7, is derived from the gauge transformation that 
connects the E and 7 gauges [l3!| and is given in Eq. ( IT7|) . Note that generically xi^'^) 
can be positive or negative. But the metric in the Einstein gauge, in either the gravity or 
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antigravity patches, has always the correct signature metric, hence, the absolute value in the 
relation for a"^ in Eq. f lT7|) . The inverse of the transformation flT7|) involves four quadrants 
in the (0^, s^) space depicted in Fig. [H and is given by 

= ± I ^V^^^^h (^) , if X > ^ ^ ^ 1 ^v^sinh (^) , if X > ^^^^ 
1 ^ Vl^sinh (^^) , if < ' [ ^^cosh , if ;^ < ■ 

As indicated in Fig. [H in the gravity sector (left/right quadrants) x{^^) is positive, while 
it is negative in the antigravity sector (top/bottom quadrants). 




FIG. 1: The gravity regime {cp'^^ > s^) is in left and right quadrants and the antigravity regime 
((/>^ < s^) is in top and bottom quadrants. 



The transformations (fT7|) and can be used by starting directly from the original 
action (|2]) without ever mentioning the Weyl symmetry or the Einstein gauge. Further- 
more, it can be used for any spacetime dependence of the fields a (x) ,aE (x), (p^ (x) , (x) 
to rewrite all equations in terms of 0^ (x) , (x) rather than a (x) ,aE (x) . If used in that 
sense then it can be considered to be an analog of what the Kruskal-Szekeres coordinates 
are to the Schwarzchild coordinates in the description of the black hole spacetime. Namely, 
they demonstrate that the spacetime regions across the horizon are smoothly geodesically 
connected. The analog here is that the field space regions (0^ ~ ■^7) > {'P'^ ~ ■^7) < 
are naturally connected at 0^ = s^; hence, the patches of spacetime x'^ in which each 
inequality is satisfied are connected geodesically, as seen in our purely time dependent solu- 
tions. Note that, in the 0^ (r) , s-y (r) form, the Friedmann equations fl28H30l) do not display 
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any singularities. Further thought suggests that the connection 0^ (x^) = (x) (which cor- 
responds to the dashed 45 degree hues in field space in Fig. [1]) can occur only at spacetimes 

at which the curvature in the Einstein frame diverges R{gE (x)) = oo. 

Inserting the above field redefinition into the standard Friedmann equations and equation 
of motion for a in the isotropic limit (derived from the action in ([2]) or f ll4p ), 
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where dot denotes derivative with respect to the conformal time r, we obtain 








+ 4:csl + Ks. 



1 



60^ + -Kcj): 



+ Pr- 



(25) 
(26) 
(27) 

(28) 
(29) 

(30) 



These are precisely the equations of motion derived from the isotropic limit of the action 



([3]) in the = 1 gauge [13|. The fields 0^, are decoupled except for the zero energy 
condition in Eq. fl30|) . These equations are valid not only when x > but also when 
X < 0, hence, they provide a smooth continuation from the gravity sector to the antigravity 
sector. The question is whether this continuation is required by the dynamics as part of 
the evolution of the universe in a geodesically complete geometry. We demonstrate that the 
generic geodesically complete solutions do require this continuation. 

Integrating the first two equations, and using the energy constraint, leads to the following 
decoupled first order differential equations 
1 > 



1 K 

2 7 ^ '"^7 ^ 2 



(31) 



with a relation between the two integration constants {E^, Eg) that reduces the unknowns 
to the single energy parameter E 



Eg = E, E^ = E + Pr- 



(32) 



Equations (13T|) are analogous to the equations satisfied by two non-relativistic particles 
with "position" coordinates, 0^ and Sj, moving independently from each other as controlled 
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by the potentials V {(j)^) = —h(l)'^ + ^<fi^ and V (s^) = cs^ + ys^, at energy levels E^j, = E + pr 
and Eg = E respectively. So, the nature of the solution and the corresponding physics is 
easily obtained through this analogy. It is sufficient to draw a picture of the potentials 
V (0^) , V (s^) and indicate the energy levels E^ = E + pr and Es = E on these pictures, 
and then let each particle begin its motion at some arbitrary points. The reader is invited 
to examine the figures in the Appendix to follow our arguments below. 

We can choose to begin the motion at tq with initial values 0^ (tq) , s (tq) that insure the 
gauge invariant factor (1 — is positive in all gauges at the initial time tq . Due to the 

time translation symmetry one of these initial values may be fixed once and for all without 
losing generality. For example, we may begin the motion somewhere on the horizontal line 
in Fig. [T] which is in the gravity sector in any gauge 

(1 - (To) (ro)) = (1 - 4 (To) (ro)) = 1 > 0. (33) 

This motion begins in the right or left quadrants and can be described initially by choosing 
the Einstein gauge (E]). The ensuing motion gives the time dependence of the generic solution 
(f)^ (r) and (r) . The solution is controlled by 6 parameters, namely the four parameters 
that define the model (6, c, K, pr) , one initial value energy parameter E and one initial value 
(f)^ (to) . We need to analyze various regions of the 6-parameter space because the motion 
can be qualitatively different in different ranges of the parameters. This is easily seen by 
staring at the pictures of the potentials in the Appendix (see [l^ for the discussion). 

The generic motion is oscillatory with each particle moving independently with indepen- 
dent oscillation periods. Each particle may pass through zero at various times independently 
from each other. Hence (0^ (r) — (r)) keeps changing sign in the 7-gauge for the generic 
solution. This shows that generically the gauge invariant factor (l — (r) (r)) , which is 
the same in every gauge (l — (r) (r)) = (1 — (r) (r)), changes sign periodically 
at times r = r„ in every gauge, where r„ is defined by the zeros of the gauge invariant factor 
computed in the 7-gauge (l — (r„) (^n)) = 0. At precisely these times the scale factor 
in the Einstein gauge vanishes as seen from Eq. f|T7|) . a"^ (r„) = ^ (0^ (r„) — (Tn)) = 0, 
and hence, this is when there is a big bang or a big crunch. In the Einstein gauge the 
generic motion must be terminated artificially at these instants of time since a|; is positive 
by definition. However, in the 0^, space the motion continues smoothly to the antigravity 
regime where 0^ (r) < (r) , which shows that the Einstein frame is geodesically incom- 
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plete for the generic motion. There exists a special subset of solutions that is geodesically 
complete in the Einstein frame without wandering into the antigravity sector, but for now 
let us continue with the generic motion. 

We have seen that even though the motion began in the gravity sector in Fig.[T]the particle 
reaches some point on the "light cone" in (0, s) space where cfp' (ri) = (ri) at r = ri in 
any gauge; it then moves smoothly into the antigravity region where 0^ (r) < (r) for some 
period of time ti < t < T2. It then turns around and passes through some other point on 
the lightcone at 0^ {T2) = {T2) at t = T2, thus reaching the gravity region again. The 
generic motion continues periodically in this way, oscillating back and forth between the 
gravity and the antigravity regions. This generic motion is represented by a curve in the 
(0^, s^) plane. Since we have the analytic solutions, we can construct the curve explicitly 
as a parametric plot 0^ (r) , (r) as shown in Fig. [2j The precise curve of the parametric 
plot is determined by the values of the 6 parameters {b, c, K, pr, E , cj) (tq)) . Generically the 
curve winds around the (0^, s^) plane since 0^ (r) , (r) are each periodic, although their 
periods are generically incommensurate. The curve becomes a closed curve if the ratio of 
the periods is a rational number. This indicates that the generic solution has repeated big 
bangs and crunches and can be cyclic when the periods are commensurate. Within each 
cycle there is a period of antigravity sandwiched between every crunch and the following big 
bang. 

The mathematical expressions of the solutions are given in the Appendix. The methods 

nn 

we used here follow those of |13l . 1141 ]. The solutions are parametrized in terms of Jacobi 
elliptic functions sn{z\m), cn{z\m), dn{z\m) [2^] in combinations chosen appropriately in 
all the relevant regions of the parameter space. For example, in a given region of parameter 
space the solution looks like 

with a similar expression for (r) , where the factors A, T, m are determined in terms of the 
6 parameters (6, c, i^', p^? E^ (tq)). The interested reader will find these expressions in our 



recent paper 



14l | for the case of = 0. The only modification to generalize to p,. 7^ is that 



previously we had used = = E, while presently we have E^ = E and E^ = E + p^.. 
Since E^ > Eg due to Pr > 0, there are more cases to investigate depending on the regions 
of the parameter space. 
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FIG. 2: The generic isotropic solution crosses back and forth through the gravity (left and right 
quadrants) and antigravity (top and bottom) quadrants. 

IV. GEODESICALLY COMPLETE BOUNCES WITHOUT ANTIGRAVITY 

Are there solutions that avoid the antigravity period in the cycle? Yes, there is a subset 
of the parameter space for which the universe completely avoids antigravity. Although 
this behavior is not generic, such solutions can be characterized as the only ones that are 
geodesically complete in only the gravity regime of the Einstein frame, which means they 
are fully described by the action ([2]) with the additional condition of geodesic completeness. 
This special subset is obtained by demanding that |0^(r)| > |s^(r)| at all times. So for 
the solutions of cf)^ (r) that oscillate between positive and negative values, at the points in 
time when (r)| vanishes (r)| must also vanish. This is possible only if the period of 
is an integer multiple of the period of s. Since there is a time translation symmetry in the 
differential equations, without losing generality we can choose that the first instance they 
both vanish is at r = 0. Hence, for 0^ (r) , (r) we must require 0^ (0) = (0) = which 
synchronizes their initial values to be both zero. This point in time is the big bang since 
we compute from Eq. (fT7|) that at this time the scale factor vanishes a^; (0) = 0. There are 
regions of parameter space that yield such solutions, but as compared to the full 6-parameter 
space it may be considered a set of measure zero. In the case of no radiation, = 0, this 
geodesically complete subset of solutions is given analytically in jl^ . 

A special example of such a solution is given in the parametric plot in Fig.[3]borrowed from 
This is a solution in a region of the parameter space where there is no radiation = 0; 
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FIG. 3: Non-generic, zero-size bounce cyclic solution that never crosses into the antigravity region. 
This figure is for 6 < and c > 0. If 6, c > the figure extends to oo in the (p direction (see 14i|). 



no curvature K = 0; special initial conditions 0^ (0) = (0) = 0; and a quantized relation 
b = —cjn^ for integer n. This quantization arises from asking the relative quantization of 
the periods for 0, s. Besides these restrictions the parameters are free to be in the regions 
c> , E > and n > 2. 

In Fig. [3], with n = 6, the fields 0^, start out initially both vanishing 0^ (0) = (0) = 
at the big bang (the arrow at the origin of the figure); then while 0^ (r) keeps growing, the 
field (r) oscillates several times until 0^ (r) reaches its maximum and turns around; then 
0^ (r) decreases to zero while (r) oscillates several times and vanishes at the same time 
as 0^. This point in time represents a big crunch. Then the motion continues smoothly to 
negative values of 0^ and repeats the same behavior of big bang then turnaround and big 
crunch. The full cycle is repeated again and again periodically which is described by Fig. 5 
in the Appendix. Note the 5 nodes in this figure are determined by the choice of the integer 
n = 6. 

When > 0, the quantization requirement for the periods puts a less severe restriction 
on the parameters (6, c, p^, E, (tq)) . Although the synchronization of initial conditions 
(0) = s (0) = and the relative quantization of the periods are still necessary, these 
conditions no longer require that b/c is quantized by itself because the additional parameter 
Pr also enters in the quantization of the periods. Instead, the parameters {b,c, K, pr, E) 
collectively are subject to one quantization condition; e.g.; the integration parameter E 
may be quantized in terms of the other four parameters plus an integer. An example of 
such a solution with radiation (but with K = for illustration) is given in the first line of 
Table I, by the parameters that satisfy ^(-^+^'') = with n = integer. This is solved by 
a quantized integration parameter En = —^r^- On the other hand when pr vanishes we 
have, ''^^^^^ = where the parameter E drops out and there is a solution like the one 
in Fig. 131 only if the parameters of the model are quantized b = —cjn^. The inclusion of 
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radiation changes the parametric plot above in a simple way: the trajectory extends further 
out in the 0-y direction as pr increases due to the higher energy in the 0^ field. 

We list below all the cases of parameter subspaces that permit purely gravity (i.e. no 
antigravity regime) geodesically complete solutions and point out the corresponding figures 
and formulas shown in the Appendix. All of these describe a universe that always remains 
in the gravity regime of the Einstein frame, and either: (i) bounces at zero size for K — 
(ii) bounces at zero size for X 7^ 0; or (iii) bounces at finite size for X > 0. These are found 
by setting 

0(O)-s(O) = O, (35) 

(which implies 5 — and then replacing — E + pr and Eg — E, instead of E^ — Eg — E, 
in the quantization of the periods. These necessary conditions cannot be satisfied for all 
the solutions given in the figures in the Appendix; the cases that are compatible with these 
conditions are indicated on the right side of Tables I,IIa,IIb,III. 

• li K — 0, there are two regimes of parameter space in which there can be a singular 
bounce without violating the null energy condition: 
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If > 0, there exist two categories of cyclic solutions, the ones that bounce at zero 
size without violating the null energy condition, and the ones that bounce at finite 
size. 

The conditions on the parameters for bouncing at finite size are 
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The conditions on the parameters for bouncing at zero size are given below. In these 
expressions /C (m) = EUipticK (m) is a well known special function that corresponds 
io the quarter period of the Jacobi elliptic functions, such as sn{z\m), with label m 
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• li K < 0, the conditions on the parameters for bouncing at zero size are 
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• In addition, for any values of 6, c, there is the special solution in which (r) = for 
all T (sitting at the s = extremum of K = \Ks^ + cs^) while (r) performs any 
motion at energy i?^ = p^. 



V. GEODESICALLY COMPLETE BOUNCES WITH ANTIGRAVITY 



In addition to the solutions described in the previous section, there are ones that are 



geodesically complete provided an antigravity regime is included (see Fig. 2). In Refs. 



17l |. we show that, when anisotropy is added to the curvature and radiation, there is a 



strong attractor behavior such that almost all solutions pass through the origin and undergo 
a period of antigravity (a loog) between each big crunch and big bang. This is illustrated in 



Fig. 4. As discussed in 16|, ll7|, the zero-size bounce solutions that evolved from crunch to 



bang in the purely gravity region in the absence of anisotropy (as listed in Tables I, lib, III) 
and illustrated in Fig. 3, as well as the other generic solutions in the Appendix illustrated in 
Fig. 2, are strongly modified near the singularity by the anisotropy, such that the trajectory 
cannot avoid the antigravity region. Furthermore, given some initial conditions, the global 
behavior of a trajectory, far away from the singularity, can also be altered even by a small 



amount of anisotropy 



171]. The finite-size bounce solutions could avoid the antigravity region 



despite anisotropy, but this may occur only in a very narrow region of parameter space. 

We are, therefore, faced with trying to understand physical phenomena in the antigravity 
regime. Since physical intuition for gravity is developed mainly in the Einstein frame, we 
begin with the Einstein gauge. When {cfp — s^) is negative, it is again possible to use the 
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FIG. 4: Comparison of a solution without anisotropy (green dotted path) and with anisotropy 
added (red thick sohd curve). An attractor mechanism caused by the anisotropy distorts the path 
so that it passes through the origin at the crunch and undergoes a loop in the antigravity region, 
through the origin again, and then re-emerging in the gravity regime. 

Weyl symmetry to choose an Einstein gauge, (p'^ — s'^ = — 1/2k^, but this is in a new domain 
of field space, namely in the top and bottom quadrants of the (0, s) space, as shown in Fig. [H 
The new (PeiSe are given by interchanging the sinh and cosh in Eq. ([6]), namely (x) = 
±^ sinh(K(T (x) /Ve), se {x) = ±^ cosh(Ka (x) /Vq), such that 0| - s| = -1/2k'^. Then 
the gauge fixed form of the action ([3]) looks like the action in Eq. ([2]), except that the first 
two terms change sign. The potential V {a) does not change sign, but it is a new function 
V (cr) , which is related to Eq. ([1]) by interchanging sinh and cosh. The metric g'^ in this 
gauge has no signature change. Hence, for matter fields, including radiation, the signs of 
their kinetic terms remain the same in the gravity and antigravity sectors. The gauge-fixed 
action in the antigravity regime looks as follows 

S - fd^x/^ah^'' ^'"^ + 1^-5.-5.- -Via)] 

J I +radiation + matter I 

Because of the sign change in the first two terms, a now looks like a ghost while the 
degree of freedom is no longer a ghost. The zero total energy constraint (Goo = ^oo Einstein 
equation) compensates for one ghost, as it did in the usual gravity regime, so there are no 
unitarity. concerns regarding the a degree of freedom. However, the other fluctuations of 
the metric, namely the spin-2 gravitons, now have the wrong sign kinetic terms. Note that 
some of the spin-2 degrees of freedom are in the form of the anisotropy fields; including them 



does not seem to show any particular instability or other unusual behavior 
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17|. 



The discussion of the relativistic harmonic oscillator, as treated in [29|, is a good toy 
model to understand how to correctly quantize the theory while maintaining unitarity when 
some degrees of freedom have the wrong sign kinetic energy. The basic technique is to 
interchange the roles of creation-annihilation operators when the kinetic en ergy has the 
wrong sign; then the resulting Fock space has only positive norms. As seen in 29| a similar 
approach also occurs in the construction of unitary representations of non-compact groups 
by using oscillators. Similarly, in the antigravity regime, the theory should be quantized 
without negative norm ghosts by interchanging the roles of creation-annihilation operators 
for gravitons. The price is that the energy of the spin-2 gravitons is unbounded below, 
so potentially there is an instability. At the linearized level, which defines a perturbative 
Hilbert space, there is no consequence. But when interactions are included, due to the 
availability of negative energy states, it may be possible to emit abundantly spin-2 gravitons 
with negative energy; however this must be accompanied with the emission of positive energy 
matter to maintain the zero energy constraint. So the theory must react in some interesting 
ways through the interactions as soon as the antigravity regime is reached. 

In the cases where there is anisotropy 16| 17|], we find that the trajectory of the antigravity 
depends on the radiation density such that, if pr increases due to the spontaneous production 
of negative energy gravitons (as noted above), the effect is to decrease the duration of the 
antigravity period and return the universe more rapidly to the big bang and a period of 
pure gravity expansion. As noted in 16|] 30|] this is an indication that the dynamics tries to 
minimize the effects of antigravity, but the details of how this works is currently cloudy. 

Of course, quantum gravity effects need to be also included. Therefore, it would be very 
interesting to study similar circumstances in the framework of string theory. To formulate 
the antigravity aspects in string theory we could use the field transformations given in 
Eqs.f ll7|l2^ . but even better would be the inclusion of the analog of the Weyl symmetry in 
the framework of string theory. 

It is worth mentioning that it seems possible to connect the state of the universe before 
the crunch to the state of the universe after the crunch by solving our classical equations 
analytically along a path in the complex r plane, such that the path completely avoids 
the antigravity regime, and also stays sufficiently far away from the singularities, so that 
quantum corrections become negligible. Such an approach is very desirable for the cyclic 
universe scenario. We will report on this type of solution in a separate paper. 
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VI. SUMMARY 



In this paper we have used analytic solutions of cosmological equations to discuss geodesic 
completeness through the big bang singularity. In the context of the path integral, our 
complete set of classical solutions provide a semi-classical approximation to the quantum 
theory. 

The computations presented in this paper mostly ignored anisotropy and used a special 
potential energy V (a) to obtain all the analytic solutions of the Friedmann equations, in a 
model that includes radiation and spatial curvature. The solutions are characterized by six 
parameters that include initial values and model parameters. We learned that the generic 
solution, in which none of the six parameters are restricted, shows that the trajectory of the 
universe goes smoothly through the crunch/bang singularities while traversing from gravity 
to antigravity spacetime patches, and doing this repeatedly in a periodic manner. The 
generic trajectory can cross the "lightcone" in field space, shown in Fig. [H at any place. 
The crossing points on the "lightcone" depend on the values of the six parameters. Although 
our general exact results are obtained in a specific model, the presence of antigravity is likely 
to occur generically in any model that is geodesically complete. Therefore, the phenomenon 
of antigravity should be considered seriously in discussing cosmology. 

We found that it is possible to avoid antigravity and still have a geodesically complete 
geometry within a smaller (but still infinite) subset of solutions (Tables I, lib, III). These are 
the only geodesically complete solutions contained totally within the traditional Einstein 
frame. One group of trajectories passes through the center of the "lightcone" repeatedly, 
resulting in a cyclic universe. These solutions, which do not violate the null energy condition, 
provide a set of examples that bouncing at zero size is possible classically in cosmological 
scenarios with or without spatial curvature. 

It should be emphasized that our new results transcend the specific simple model above. 
The phenomena we have found should also be expected generically in supergravity theories 
coupled to matter whose formulation include a similar factor that multiplies R (g) . In su- 
pergravity, that factor is related to the Kahler potential, and this factor, combined with the 
usual Einstein-Hilbert term, is not generally positive definite 3l|. In fact, in a gauge that 
we call the supergravity gauge, or c-gauge, in which (j){x) is set to a constant (j)o 3], 
our term (0q — s'^)R{g) reduces precisely to the familiar form in supergravity including a 
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Kahler-like potential. In the past, it was assumed that the overall factor is positive and in- 
vestigations of supergravity proceeded only in the positive regime. A discussion of the field 
space in the positive sector for general N'=2 supergravity can be found in [32]. However, our 
results suggest that generically the overall factor can and will change sign dynamically, in 
every gauge, and therefore antigravity sectors similar to our discussion in this paper should 
be expected in typical supergravity theories. This is illustrated with an example in jl^. 




Until better understood in the context of quantum gravity, or string theory, our results 
should be considered to be a first pass for the types of new questions they raise and the 
answers they provide. 

Much remains to be understood, including quantum gravity and string theory effects, 
but it is clear that previously unsuspected phenomena, including antigravity, come into play 
classically close to the cosmological singularity. The technical tools to study such issues 
in the context of a full quantum theory of gravity are yet to be developed. This is an 
important challenge to the theory community, since the results have profound implications 
for both fundamental physics and our understanding of the origin, evolution and future of 
the universe. 
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Appendix A: The Analytic Solutions 

The intuitive approach for solving the Friedmann equations in the separable form 
Eqs.(l28]— I 301) was described in Sec. IIIII In Eqs.( l30|32l) we showed that a first integral of 
these second order equations takes the form of the analog problem of a particle in a poten- 
tial, separately for the "particles" 0^, s^. Namely, \(j?^ + V^ = E^, and \s^^ + Vs = Es, with 
the potentials = —bcj)^ + and K = cs^ + with the energy constraint 
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generally solved by Eg = E and Ef/, = E + pr for any E. In this appendix we list the complete 
set of analytic solutions to these equations in all possible regions of the six parameter space 
(6, c, K, pr, E, (f) (to)) . In this Appendix, we restrict ourselves to the case 6 + c > 0, which 
corresponds in the Einstein frame to a potential V (cr) in Eq. ([T]) that is bounded from below. 
The generalization to unstable potentials is straightforward. 

Figs. 5-29 below represent the potentials (solid line) and Vs (dashed line) and the 
energy levels E^, Eg] they are drawn in the various regions of parameter space. The 0- level 
E^ = E + pr is higher than the s-level Eg = E because the radiation energy density pr is 
positive. The parameter E is allowed to slide vertically within the physical range permitted 
for the s-particle's motion in the potential Vg. From the figures alone one can obtain the 
intuitive solution by invoking the analogy of a particle moving in a potential for either or 
(p-y. Next to each figure we give the corresponding analytic solution to Eq.fl28]- 1 30p . We did 
not include separately some trivial cases such as the case when the fields sit still at the top 
or bottom of the potential extrema, or in the cases of 6 = or c = where Eqs.(!28]— I 30p 
become linear differential equations with simple trigonometric solutions. These cases are 
recovered in the appropriate limit of the expressions below. 

In the analytic expressions below instead of the free parameter (p (tq) for an arbitrary 
initial value, we have inserted the arbitrary phase shift parameter S. Note also that we 
have used the time translation symmetry of the equations to choose a specific value for the 
initial value for s-y (r) at r = 0. In this way we give here the generic solutions that describe 
all possible geodesically complete trajectories of the universe, including those that move 
between the gravity and antigravity regions. 

A subset of solutions stay only in the gravity sector with geodesically complete trajec- 
tories. They are obtained by putting constraints on the parameters. These include setting 
6 = (synchronized initial conditions for 0, s) and requiring a quantization of the periods 
of relative to the period of s. The corresponding parameter space is given in section ( II Vp . 

The reader can verify that the following expressions solve the differential equations and 
that the plots of (r) , s (r) as functions of r correspond to the motion intuitively expected 
for a particle in the corresponding potential at the given energy level. To verify the solution 
the following properties of the Jacobi elliptic functions sn {z\m) , cn {z\m) and dn {z\m) are 



useful 



28| . The derivative of Jacobi elliptic functions are given in terms of expressions 
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somewhat similar to those for trigonometric functions 

—sn {z\m) — cn {z\m) x dn {z\m) , (Al) 
dz 

-^cn {z\m) = —sn {z\m) x dn {z\m) , (A2) 
dz 

-j-dn (zlm) — —m x sn {z\m) x cn {z\m) . (A3) 
dz 

They also satisfy quadratic relations, such as 

{sn {z\m))'^ + {cn {z\m))'^ = 1] m{sn{z\m))'^ + {dn{z\m))'^ = 1. (A4) 

When m is a real number in the range — oo < m < 1, the function sn {z\m) oscillates 
between the values —1 to +1, similar to the trigonometric function sin{z), vanishing at 
z — )■ 0, and reaching a maximum at the quarter period z — K {m), where K {m) is the 
elliptic integral as a function of m. Prom the quadratic relations above, it is deduced that 
the behavior of cn{z\m) is that it oscillates similar to a cosine cos (2;) , while dn{z\m) 
oscillates between the positive values +1 and and (1 — nif'^'^ . When m > 1 the behavior is 
still oscillatory but quite different than sin (2;) , cos (z) ,etc.. However, it is possible to use 
identities to rewrite the solution in terms of sn {z\m') , cn {z\m') , dn {z\m') where m' = 1 — m 
is again in the range — 00 < m' < 1. We have used such identities so that all the m values 
that appear in our solutions below are in the range — 00 < m < 1. Then the reader can get a 
feehng of the behavior of the solutions by the analogy to trigonometric sin {z) , cos {z) ,etc.. 
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We begin with the K = cases; there are five different regions for the remaining param- 
eters as hsted in Figs. (5-9) 
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Note that as the parameters b, c, E, {E + pr) change signs the corresponding solutions 
and physical behaviors change qualitatively. Nevertheless, the mathematical expressions in 
Figs. (5-9) are related to each other by the following rules for analj^ic continuation, where x 
is real 



So, it is possible to write a single formula to cover all the solutions (such as the formulas in 
Fig. 9, modified with appropriate absolute signs) 



and then analytically continue the argument of the Jacobi elliptic functions to obtain the 
other expressions. In this form all signs of the parameters 6, c, E, {E + pr) are permitted, 
thus capturing the physical behavior of all corresponding regions of parameter space with 
a single expression (when K — 0). Under these flips of signs the functions 0, s remain real 
even though the argument of the function is complex. This unified version is convenient to 
feed it to a computer to obtain plots of the solutions. 

There is a similar analytic continuation for the cases with nonzero curvature given below, 
but the formulas for analytic continuation are considerably more involved, so we will not 
bother to discuss them. 




(A5) 



(A6) 




(A7) 
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Next we have the K > cases, with eleven combinations which are hsted in Figs. (10-19) 
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FIG. 14 

6 > 0, c < 
Es<0<E^< 
10(0)1 > 



46' 



166 

.(0)1 > 



ibT'i 



cnl ^^\m^ 



m 



m^ = i(l-irT|) <0 
T^ = {K^-16b{E + pr)y'^' 
\{l-KT^)<\ 



4|c|T| cn(^|m,)' 



T, = {K^ + lQ\c\\E\) 



-1/4 




FIG. 15 

6 > 0, c < 0, 

.ft' 



l</>(0)|< 1^(0)1 >V^ 



iiT^-l 
2bT^ 



sn 



KTl - 1< 1 

K 



+ 



4 



- 46 (E + p^) 



4|c|T| cn(^|r 



i (1 - KT^) < 



Ts = {K^ - 16 |c| E) 



1 



-1/4 




FIG. 16 

6 > 0, c < 0, E^>Es 



E^ > 



El 

166' 



E., < 



16|c| 



^ ^/4|c| 



0(r 
s (r 



' = (64& {E + p,))"'^' 
i(l-KTf)<0 



4|c|r| cn(^|m,)^ 



= {K^ - 16 |c| E) 



-1/4 




FIG. 17 

6 > 0, c < 0, E^>Es 

^0 > ik' < < 4^ 



k(o)l< J4 



l-J m,^ = 1 + iTTl < 1 



' = (646 (E + p.))"'/^ 
isTT? - 1< 1 



2 ^ 



f -4|c|E 



-1/2 



33 



_E(P_ 



FIG. 18 

b>0, c < 0, 
Es> Es> 0, Es < 



El 

16b' 



10(0)1 > (0)1 <Ji 



cn 



^1 



'''' ^ 



m. 



= 


= 1(1- 




< 


T^ = 




m{E 




ITLs = 




- 1 < 1 










\ -1/2 


T = 

s 




\J-- 


4|c|E 1 



-^s 

/ 

FIG. 19 

6 > 0, c < 0, 
E(i, > Eg, < Etf) < 



El 

166' 



10(0)1 < Vi' 1^(0)1 <a/^ 



/T^Tj-l 



= KTl-\<\ 



KTl-\ fj_\ 
2\c\Ts \Ts 1"^' 



= fsTT^ - 1 < 1 

-1/2 


















' ^+cn[^\m,) ' = (646 {E + pr))-'" 


/ / w 


FIG. 20 




6 > 0, c < 0, 




E^ > Es, Es > ^Ij^ 





34 



Finally, for K < there are nine combinations which are listed in Figs. (21-29) 
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